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Abstract 

In this paper, we revisit the enumeration of directed animals using gas models. We show that 
there exists a natural construction of random directed animals on any directed graph together 
with a particle system that explains at the level of objects the formal link known between the 
density of the gas model and the generating function of directed animals counted according 
to the area. This provides some new methods to compute the generating function of directed 
animals counted according to area, and leads in the particular case of the square lattice to 
new combinatorial results and questions. A model of gas related to directed animals counted 
according to area and perimeter on any directed graph is also exhibited. 



^ Introduction 

> : 

^ . . 

O^.l Directed animals on a directed graph 

' Let G = {V, E) be a connected graph with set of vertices V and set of edges E. An animal A 
^vQn G is a subset of V such that between any two vertices u and v m. there is a path in G having 
■~-«»all its vertices in A. The vertices of A are called cells and the number of cells is denoted |j4| and 
-galled the area of ^. A neighbor n of A is a vertex of G which is not in A and such that there 
^xists e £ E between u and a vertex v{u) of A. The perimeter of A, denoted by V{A), is the set of 
j^eighbors of A and its cardinality is denoted by 1^(^)1. 

Directed animals (DA) are animals built on a directed graph G : 

j^efinition 1.1 Let A and S be two subsets of V , with finite or infinite cardinalities. We say 
that A is a DA with source S, if S is a subset of A such that any vertex of A can be reached from 
an element of S through a directed path having all its vertices in A. 

In the setting of DA, the definition of cells and area are the same as in the case of animals, but 
the notion of neighbor is changed since the edge {v{u),u) is required to be a directed edge of G. If 
there is a directed edge from v to u m G, then u is said to be a child of v, and v to be the father 
of u; this induces a notion of descendant and ancestor. Each node of V{A) has at least one father 
in A. In this paper, we deal only with DA built on graphs having some suitable properties : 

Definition 1.2 A directed graph G = {V,E) is said to be agreeable if 

(A) G does not contain multiple edges, 

(B) the graph G has no directed cycles, 

(C) the number of children of each node is finite. 



Condition (A) is needed to avoid some useless discussions in the Theorems. Condition (C) is needed 
to have a finite number of DA with a given area, for all sources. Notice that an agreeable graph 
is not necessarily connected, and is not necessarily locally finite (some nodes may have an infinite 
incoming degree). Even if never recalled in the statements of the results, V is supposed to be finite 
or countable. 

As examples, (finite or infinite) trees and forests are agreeable graphs, the square lattice Sq = 
directed in such a way that the vertex (x, y) has as children (x, y + 1) and (x + 1, y + 1) is agreeable 
(as well as all usual directed lattices). 

A subset S" of y is said to be free if for any x,y £ S, x ^ y, x is not an ancestor of y. For any DA 
A, the set 

S{A) := {x, X £ A,x has no father in A} 

is a free subset of V and is the unique minimal source of A according to the inclusion partial order 
(it is also the intersection of all possibles sources of A). 




Figure 1: On the first line, three examples of agreeable graphs on which the directed edges are directed 
upwards. The first example is a tree, the second is the square lattice and the third a "non-layered" 
agreeable graph. On the second line are represented some DA on these graphs; filled points are the 
cells, the crosses are perimeter sites, and the surrounded points are the minimal sources of these DA. 



We denote by A{S,G) (or more simply A{S) when no confusion on G is possible), the set of 
finite or infinite DA on G with source S, and by (or G^) the generating function (GF) of finite 
DA counted according to the area : 

Gf(x) :=xl^l + ... 

where . . . stands for a sum of monomials whose degree are at least \S\ + 1. We will need also 
sometimes to consider DA having their sources included (or equal) in a given set S. For such a DA, 
S is called an over-source . We denote by A{S, G) (or A{S)) the set of DA having S as over-source, 
and by Gg (or Gs) its GF {G^{x) = 1 + |5|x + ...). Finally, we set Gg = G^ = 1. 



When dealing with elements A of A{S, G), the set are also considered as (special) perimeter 
sites, and we set Vs{A) =V{A)yj {S\A). 

We introduce a notion of particle systems, or gas occupation, on a graph : 

Definition 1.3 Let G = (V, E) be a graph. A particle system on G, or gas occupation on G, is a 
map C from V to {0, 1}. The vertices v such that C{v) = 1 are said to be occupied, the others are 
said to be empty. 

In the physic literature, a hard particle model (or occupation) on a graph is a (probability model 
of) gas occupation C with the additional constraint that the occupied sites are not neighbors in G. 
When the gas is random, we call density of C in a vertex x, the quantity ¥{C{x) = 1). 

1.2 Contents 

The equivalence of enumeration of DA on lattices according to the area and solving hard par- 
ticle models has been a key point in the study of DA from its very beginning with the work of 
physicists. Hakim & Nadal [12], Nadal, Derrida, &; Vannimenus [14], Dhar [8], [10] and later, thanks 
to some generalizations due to Bousquet-Melou [3] and Bousquet-Melou & Conway [4]. The idea 
is to consider on the same lattice where are biiilt the DA an ad hoc random gas models owning 
the same "recursive decomposition" as the DA up to a slight change of variables. This recursive 
decomposition is done using a decomposition of the lattice itself by layers (see Section 5.2, where 
the historical approach used by Bousquct-Mclou [3] is detailed). Dhar [9], using statistical me- 
chanics techniques, explained that in special cases there exists an exact equivalence between the 
enumeration of DA in d dimensions and the computation of the free-energy of a (d — l)-dimensional 
lattice gas. A combinatorial explanation using heaps of pieces was proposed by Viennot [16] (see 
also Betrema Sz Penaud [6] for a pedagogical and detailed exposition). 

To solve the equation involving the GF of DA obtained by this recursive decomposition some 
properties of the gas model are used. The gas model is a stochastic process indexed by the lattice 
having in the tractable cases some nice Markovian-type properties on the layers. The rigorous 
arguments given in [3, 4] avoid the construction of the gas model on the whole lattice as done by 
Dhar : the gas model is defined on the layers of a cylinder (a lattice having some cyclical boundary 
conditions), and the transition allowing to pass from a layer to the following one are of Markovian 
type. Bousquet-Melou [3] finds an explicit solution of the gas process on a layer (in the square 
lattice case, in the triangular lattice case, and in other lattices in the joint work with Conway). 
Then, the computation of the density of the gas distribution is explicitly solved using that the 
number of configurations on such layers is finite, leading to rational GF. The GF of DA on the 
entire lattice (without the cyclical condition) is then obtained by a formal passage to the limit. 

In this paper we revisit the relation between the enumeration of DA on a lattice, and more 
generally on any agreeable graph, and the computation of the density of a model of gas. Our 
construction coincides with that used by Dhar or Bousquet-Melou in their works. In Section 2, we 
explain how the usual construction of random DA on a graph G, using a Bernoulli coloring of the 



vertices of G, allows to define in the same time a random model of gas (that we qualified to be of 
type 1). Here the construction is not done in "parallel" as in the works previously cited but on the 
same probability space; this provides a coupling of these objects. Using this coupling, we provide 
a general explanation of the fact that the GF of DA with source {x} counted according to the area 
on any agreeable graph equals the density of the associated gas at vertex x, up to a simple change 
of variables (Theorem 2.6). This explanation is not of the same nature than in the previously cited 
work: the link between the density and the GF is not only formal but is explained combinatorially 
at the level of the DA (Section 2). Moreover, the construction of the gas model can be done also 
on (non-regular) infinite graph in a rigorous manner and no passage to the limit is needed. 

In Section 3, we revisit the study of DA on the square lattice; the new description of the gas 
model on the whole lattice allows us to provide a description of the gas model on a line (Theorem 
3.3). On this line the gas model is a Markov chain which is identified. We provide then a new 
way to compute the GF of DA counted according to the area (Theorem 3.3). This extends to the 
enumeration of DA with any source on a line (Proposition 3.6); till now only the case with compact 
sources was known. We explain also how to compute the GF of DA with sources that are not 
contained in a line (Remark 3.8) and provide an example (Proposition 3.7). 

In Section 4, we present an other model of gas, that we qualify to be of type 2. The density of 
this gas model is related to the GF of DA counted according to the area and perimeter (Theorem 
4.3). This construction explained once again at the level of object on any agreeable graph a relation 
used by Bousquet-Melou [3] in a formal way on the square lattice. Even if we haven't find any 
deep application to this construction, we think that it provides an interesting generic approach to 
the computation of the GF of DA according to the area and perimeter, and it should lead to new 
results in the future. 

Some other references concerning DA on lattices 

One finds in the literature numerous works concerning the enumeration of DA on lattices, most 
of them avoids gas model considerations. We don't want to be exhaustive here (we send the reader 
to Bousquct-Melou [3], Viennot [15, 16] and references therein), but we would like to indicate some 
combinatorial works directly related to this paper. It is interesting to notice that an important part 
of the papers cited below are combinatorial proofs of results found before using gas techniques. 

First, we refer to Viennot [16] (and Betrema &; Penaud [6]) for the algebraico-combinatorial 
relation between DA and heaps of pieces. This powerful point of view having some applications 
everywhere in the combinatorics, allows to compute the GF of DA on the triangular lattice, and 
by a change of variables on the square lattice. A direct combinatorial enumeration of DA on the 
square lattice has been done by Betrema &; Penaud [5]; they found a bijection with a family of 
trees :" les arbres guingois" . 

Heap of pieces techniques have been used by Corteel, Denise &: Gouyou-Beauchamps [7] to give 
a combinatorial enumeration of DA on some lattices, first counted by Bousquet-Melou & Conway 
[4] using gas model (of type 1). Viennot and Gouyou-Bcauchamps [11] provide a bijection between 
DA with compact sources on the square lattice and certain paths in the plane; they are able to 



enumerate these DA. Barcucci k, al. [2] studied DA on the square and triangular lattices with 
the help of the ECO method. They found some relations with permutations with some forbidden 
subsequences and a family of trees. 

2 Simultaneous construction of DA and gas model of type 1 

In this part, we construct on any agreeable graph G a probability space on which are well- 
defined a model of gas - that we qualified to be "of type 1" - and a notion of random DA. The 
relation between the density of the gas model of type 1 and the GF of DA counted according to 
the area will then be explained at the level of objects on this probability space. 

2.1 Construction of DA 

Let G = {V,E) be an agreeable graph. We introduce a random coloring of V by the two 
colors a and b. We need to be a little bit formal here since when G is infinite the existence 
of a probability space where such a construction is possible is not so obvious. We consider the 
probability space J7 = {a,b}^ endowed with the c-fields J- generated by the subsets having a 
finite cardinality (the cylinders in the usual probability terminology) and endowed by the measure 
product Fp = (p5a + (1 — p)Si,)^^ , where 5a is the standard Dirac measure on {a}; let C = {Cx)xeV 
be the identity map on V. We denote by Ep the expectation under Pp. 

In other words C is a random coloring of V : under ¥p the random variables Cx giving the color 
of the vertex xofV are independent and take the value a and b with probability p and 1 — p. 

Let uj e and 5 be a subset of V. We denote by S,{uj) = {x,x G S, Cx{oi)) = a} the subset 
of S of vertices with color a. We denote by A^{uj) the maximal DA for the inclusion partial order 
with source S,{u) and whose cells are the vertices x such that Cx{uj) = a that can be reached from 
S,{uj) by an a-colored path. By construction the perimeter sites of A^{uj) are 6-colored (see Fig. 
2). 




S 



Figure 2: The DA A^{uj) is the set of gray cells. 



Proposition 2.1 Let G = {V,E) be an agreeable graph. 

{i) Let S be a free subset of V and B G A{S, G) be a finite DA with source S. We have 

Pp(A^ = 5)=pl^l(l-;>)l^(^)l. 



(u) Let B G A{S, G) he a finite DA with over-source the free set S. We have 



Proof {i) First {lo) = B \i and only \iCx = a for aWx e B and \iCx = b for all x G V{B). Then 
{A'^ = i?} is ^-measurable since it depends only on a finite number of cells, and the conclusion 
follows. For (ii), we use moreover that Vs{A) := V{A) U (5 \ A). □ 

2.2 Directed animals and percolation 

When G is an infinite graph and |5| > 1, under Pp the random DA A'^ may be infinite with 
positive probability. The probability to have an infinite DA with source S is also that of the 
directed sites percolation starting from S where the cells of the percolation cluster are the vertices 
with color a reachable from S by an a-colored directed path. 

Denote by p^^^ the threshold for the existence of an infinite DA with positive probability : 



Most of the results of the present paper are valid only when p < ^J^rif "^^^ threshold p^^^^ is in 
general difficult to compute, but here is a simple sufficient condition on G for which p^^.^^ > 0. 

Proposition 2.2 Let G be a agreeable graph such that the maximum number of children of its 
vertices is bounded by K. Then for any finite subset S of G, p^j-u — '^/K. 

A proof of that result is given in the Appendix. Also in the Appendix, Comment 5.1 provides a 
graph in which pcrit = 0. 

We recall two results giving some insight on the percolation probabilities (and easy to prove). 
• Let Si and S2 be two subsets of V. We have pf^'^^'^" = min (pf^ , p^^. -\ . 



• Let Pcrit = '^'^'^vev Pcrif P ^ [Oi Pcrit)) Under Pp, almost surely (a.s.) all DA in G having 

a finite source are finite. Indeed a countable intersection of a.s. events is a.s.. 

2.3 Construction of the gas model of type 1 

The construction of the gas model of type 1, Proposition 2.4 and the Nim game construction 
presented in this section are generalizations and formalization of the work of the first author [13, 
section 1.4]. 

Let us build a gas model X on an agreeable graph G = {y,E) (see Definition 1.3). This 
construction takes place on the probability space O introduced in Section 2.1, and X is defined 
thanks to the random coloring C. 

For any x and a; G denote by 



Pcrit 



sup{p, Pp(|A^ 



< +00) = 1}. 





(1) 



= 1 



n (i-^c(a.)). 



(2) 



c: children of x 



If X has no children the product in (1) is empty, and as usual, we set its value to 1. 

For any co and any x eV, Xx{u!) is to be interpreted as the gas occupation in the vertex x. 

We have to investigate when the recursive definition giving Xx{lo) is correct, that is when it 
allows to indeed compute a value Xx{u!) (see Fig. 3). 

• When Cx{io) = b then Xx{ol>) = (there is no problem to define Xx{io)). 

• When Cx{uj) = a, to compute Xx{co), knowing all the values Xy{uj) for y child of x is sufficient; 
their values are given by the same rule. By successive iterations a "DA of calculus" growths. This 
DA is exactly A^^\uj). Indeed, on each cell of A^^^(u;) - whose color are a by construction - the 
following computation is done 

XxH^ n (l-^c(a;)); 

c: children of x 

since Cy{u) = b for any perimeter sites of A{^}(a;), Xy{Lo) = on 7'(A{^})(a;)). Then, one sees 
that Xx{lo) is well defined if A^'^^{uj) is finite because this recursive computation of Xx{u}) ends. 

Moreover, in this case, if A^^^(u;) = A is finite the value Xx{uj) is a deterministic function of 
A, that we denote by Xx(A), since the only data in this calculus is the geometry of A (the map 
Xx is defined only on finite DA with source x). The maps {Xx)xeV satisfies then for simple reasons 
the following decomposition. Let v be a vertex in G, A^""^ a finite DA with source v, and denote 
by fi, . . . , Vd the children of v in G, and A'^^''^^, . . . , A^'"''-^ be the maximal DA included in A^'"^ with 
over-source {vi}, . . . , {vd} respectively. Then 

X. (aM) = l|^i„>|>on (l - XvM'^'"'^)) ■ (3) 

i=l 

In the same vein, assume that a finite free subset S of G is given. The vector {Xx{uj))xeS giving 
the gas occupation on S is also a deterministic function of the DA A'^(a;). 

Remark 2.3 In the case where Cx{c^) = a but \A^^^{uj)\ is infinite, the computation of Xx{uj) may 
also ends within a finite number of steps (using the fact that the product Yl^. children ofxi^ ~ Ac(i^)) 
is known to be when one of its terms is null, which can be the case on a finite sub animal of 
A^{lo) ). But, in this case the value of Xx returned by this procedure is inadequate and must be not 
considered in the analysis since the a.s. finitness of A^ is crucial in the proofs. In the following 
we will restrict ourselves to p < pcrit in order to avoid infinite DA A^{u) with probability 1. 

By the previous consideration we may conclude by the following proposition. 

Proposition 2.4 Let G = {V,E) be an agreeable graph, x E V and p G [0,^^^^^). Under ¥p the 
random variable Xx is a.s. well defined by (1), and 

Ep(A,) = ¥p{Xx = 1) = Ep(x.(AW)) = pl^'(l 

AeA{{x},G),\A\<+oo,Xx{A)=l 

where A{{x},G) has been defined in Section 1.1. 



Proof. For p G [0,pi^/J, the DA A{^> is a.s. finite, and then a.s. Xx{-^^^^) = X^- The three 
equaUties are straightforward. One also checks that is indeed a random variable, that is J-- 
measurable. For i £ {0, 1}, 

X-^{{i}) = [J{A,\A\ =n, DA with source {x},Xx{A^''^) =i}. 

n>0 

This sets are measurable since they depend only on a finite number of cells. The conclusion follows 
the fact that a countable union of measurable sets is measurable. □ 




Figure 3: On the first column on top, a random coloring. Below, the family of DA derived from it. 
On the second column on top, the beginning of the computation of the gas occupation. 'V stands for 
the places where the calculus ^ He- children of ~ "'^c) must be done. Below, the gas occupation 
has been computed. 

Remark 2.5 In the case where AS^^{ijj) = A is finite, it may be convenient for the reader to see the 
computation of Xx{io) as the result of a Nim type game with two players, player(O) and player(l) 
on A according to the following rules : 

- the first player that can not play is the looser. 

- at time 0, player (0) places a token on x (if A is empty, player (0) is the looser). 

- then, the players move in turn the token upward in A. At time i, player{i mod 2) takes the 
token from where it is, say in v, and can move it in u, if u is in A and if {v, u) is a directed edge 
ofG. 

The fact is that Xxi^) = 1 iff plo,yer(0) has a strategy to win against all defense of player (1). 
We let the reader proves this property as an exercise. Indication: proceed to the computation of 
Xx{A) by the last moves of the players. 

One of the main aim of this paper is to provide an explication at the level of objects of the 
following theorem 



Theorem 2.6 Let G = {V, E) he an agreeable graph, x he a vertex of V , Rx he the radius of 
convergence of G^^y and p E [0,Rx)- We have 

EpiXx) = -Gf,}(-p). (4) 

We will see in Lemma 2.9 that in any graph Rx ^ PcrW 

This relation between the density of a gas model and the GF of DA in the case of lattice graph 
is first discovered by Dhar [9] in the square lattice case; it is then made rigorous and generalized 

by Bousquct-Mcloii [3. 4]. In each case, a formula similar to (4) is obtained, but there the equality 
is only formal : Kp{X^) and — G|'^_|(— p) are shown to be formal series satisfying the same recursive 
decomposition. We want also to point out that in [3, 4], the gas model is studied on a cylinder and 
some arguments using the finiteness of the number of states, and the convergence of some Markov 
chains are used (see Section 5.2). Here these steps are not needed since the construction of the gas 
model is done "on the right graph at once" and the construction provides also the gas under "the 
stationary regime" if the graph is a lattice. We want also to stress on the fact that Theorem 2.6 
holds on any agreeable graph and not only on lattices. 

In order to express our relation in the level of objects, we will rewrite the right hand side of (4) 
under the form of an expectation, in order to make more apparent that Kp{Xx) and ~G^|(— p) 
are both sum on weighted DA (the weight of a DA A being simply p'^{l — its probability 

on O). But first, we need to introduce the notion of sub animal. 

Definition 2.7 Let A and A' be two DA. We say that A is a suh animal (suh-DA) of A' (we write 
A -< A'), if A is included in A' and if S{A) is included in S{A'). 

For any DA A with minimal source S{A) = S, and p € (0,p^j^), we have 

A' : A^A',S{A')=S A' : A^A',S{A')=S 

Indeed, is the probability of the event {lu,Xx{u;) = a for any x & A}, and this latter equals 
{lo, a is a sub-DA of A^{lo)}. We have 

Proposition 2.8 Let G = (V, E) he an agreeahle graph. For any x E V and any p € [0, Rx), 

-Gf,}(-p)=Ep(D,(AW)) 

where for any A with source {x}, 

B : B~iA,S{B)={x} 

is the difference between the number of sub-DA of A having an odd number of cells and those having 
an ( non zero ) even number of cells. 

Before proving this proposition, we establish the following Lemma 



Lemma 2.9 For any agreeable graph G = (V, E) and x in V , 

Rx<ptl (6) 

Proof. Since for any A with source is a sub-DA of A^^^) (where A^^^ may be 

infinite), 

^pl^l = ^F{A ^ A^^>) = = Number of sub-DA of A^). (7) 

A A A 

If p G [0,Rx) then these quantities are finite. This impUes that Pp a.s. the number of sub-DA of 
A^^^ is finite, which in turn imphes that A^^^ is Pp a.s. finite, which finally yields p < p^\- □ 

Proof of Proposition 2.8. It relies on a permutation of sum symbols. For any p G [0, Rx), 

A,S{A)={x} 

Since by Lemma 2.9 we have p < plrip this is equal to 

E ( E p^^'Ki-pf^^'H-irA. (8) 

A : S{A)={x} \A' : A-iA' ,S{A')={x} / 

This sum converges absolutely when p < since with an absolute value it is the mean of the 
number of Sub DA of A^^^ (which is finite when p < R^, see the proof of Lemma 2.9). By the 
Fubini's theorem the double sum in (8) is 

^ pl^'l(l -p)l^(^')l (z),(^')) = Ep (z).(AW)) . □ 

A' ■ S{A')={x} 

Remark 2.10 • In general the inequality (6) is strict, since 

Rx = suTp{p,Ep{Number of sub-DA o/ A^^^) < +00} when Pcrit = sup{p,Pp(|A^^^| < +00) = 1}. 

In the case of the square lattice Rx = 1/3 when Pcrit is larger than 1/2 (by Proposition 2.2) and is 
expected to be around 0.6. 

• Proposition 2.8 is a relation between two quantities that are not defined in general for the same 
values of p. Viewed as formal series in p, they are equal. But viewed as real functions of p, this 
is not so simple, and this is rather interesting. This is due to the classical well known fact that a 
sum which is not absolutely converging is well defined only if an order of summation of its terms is 
provided. When p > Rx the fact that ~G^|(— p) is a non convergent sum comes from the fact that 
the order of summation is given by the sizes of the DA. The second sum 5^^p''^'(l— p)'^^"^^' {Dx{A)) 
contains the same terms as ~G^|(— p) as one may guess by the permutation of the sum symbols 
authorized for some p by the Fubini's theorem. But these terms appear in some sense by groups, 
and it is a fact that this sum converges more easily than — (— p) (as said in the first point). 



We saw that -G^^y{-p) = Ep (Z)j;(A{^>)) for p < < pl^l and the density of the gas Ep{Xx) 
is equal to Ep{xx{-^^^^)), for p < p^\- Here is the exphcation of Theorem 2.6 at the level of object : 

Theorem 2.11 For any finite DA A with source v, we have Dy{A) = Xv{-^)- 

This allows to deduce that Theorem 2.6 holds true, since a.s. when p < pj.^jf, the random DA A^''^ 
is a.s. finite, and then X^{A^'"^) = Xv{-^^"^) = Dy{A^'"^) a.s., and then these variables have the 
same expectation (notice that this also implies that Dy(A) takes its values in {0, 1}, which is not 
necessarily obvious). 

In order to prove Theorem 2.11, we will show that Dy owns the same recursive decomposition 
as Xv given in (3); this is done via the introduction of a notion of embedding of trees in DA. An 
heuristic is given in Section 2.5. 

2.4 Embedded trees 

Let G = {V, E) be an agreeable graph. The set V being at most countable we assume from now 
on that an order denoted by < is given on V. This order induces an order among the children of 
a given vertex in G. Since we are to "canonically" embed some ordered trees in G we need also to 
define a suitable ordering of the nodes of those trees; this is inspired from the Neveu's definition of 
trees. 

+ 0O 

Let N = {1, 2, 3, . . . } be the set of non-negative integers and W = {0} U N* the set of finite 

i=l 

words on the alphabet N, where denotes the empty word. We define the concatenation product 
of two words u = ui . . . Uk and v = Vi . . .vi oi W hy uv := ui . . . u^vi ■ ■ - Vi; the empty word is 
the neutral element for this operation : 0u = u0 for any u eW. 

Definition 2.12 We denote by tree a subset t of W such that € t and if u = ui . . . ^ t for 
some k > 1 then ui . . . Uk-i G t. In other word if a word u is in t, its prefixes are also in t (see 
Fig. 4). 

The set of trees is denoted by T. 

In the combinatorial literature, what is called tree here is sometimes viewed as a depth first 
traversal encoding of trees. 

The vocabulary attached to trees arc as usual: is called the root, the strict prefixes of u G t 
are called ancestors of u, and m is a descendant of its ancestors. The elements of t are called nodes, 
and for G t, \u\ stands for the number of letters in u (by convention |0| = 0) and is called the 
depth of n. If u = -ui . . . Uki then ui . . . is the father of u, and v = u\. . . Ukj G t for i ^ j, we 
say that u and v are brothers. For any u = ui . . .Uk e t, we denote by C„(t) = {ui, z G N} fl t the 
set of children of in t. The size of t denoted by |t| is the cardinality of t. 

Let t be a tree and i G C0(t) a child of the root in t. For any z G N, we denote by t* the set : 



:= {v,u := iv G t}. 



2,5 2,6 7,5 1,1 1,2 3,1 




Figure 4: The trees admit some usual representation as embedded figure in the plane, mapping the 
order between brothers into the usual order between the abscissas. In the usual Neveu's convention, 
there is an additional axiom (if v = ui . . . Ukj G t for some j > 2 then v — ui . . . Uk{j — 1) G t) which 
ensures that the "branching structure of the tree" characterizes the tree (for Neveu, only the second 
drawing represents a tree). This is not the case, here. 

It is easy to see that t* is a tree, and that it is obtained from t by deleting all the nodes that are 
not descendant of and "rerooted" va. i : i and the descendants of i are kept but they lose they 
first letter in order to form a tree (see Fig. 5). The tree t* is called fringe subtree in the literature. 




Figure 5: A representation of the tree t = {0, 1, 3, 4, 6, (1, 1), (1, 4), (4, 3), (6, 2), (1, 4, 3), (4, 3, 1), (4, 3, 3)}. 
The second tree is t"*. We have C2i(t) = {1,3, 4, 6}. On the third figure is drawn the embedding 7r^(t) 
on an agreeable graph, where the edges are directed upward. 

Definition 2.13 Let G = {V, E) he an agreeable graph, v £ V and t a tree. Let vi, . . . ,Vd be the 

d = dv children of v in G sorted according to < . We say that t is embeddable in G at v ift = {0} 

or if for any i G C^it), V is embeddable at vi. If t is embeddable in G at v, we denote by 7r^(t) its 
embedding in G at v defined by : 

7r„(t) = Hu U ^-.(*')- 

ieC0 (t) 

An illustration of this embedding is given in Fig. 5. 

In other words, 7r^(t) is the subset of G obtained by drawing t on G according to the following 
rules: first draw the root of t on f (in other word 7rt,({0}) = v). Each branch of t is a succession 
of nodes = ui . . . where is a prefix of a^+i (and oq = 0). Then draw a^+i in G in such a 
way that ak+i is the Uk+i child of (the edge (7r„(aA:), vr„(afc+i)) of G is the Uk+ith. edge starting 
from 7r„(afc)) (see Fig. 5). Hence each branch of t is finally embedded in a simple path of G issued 
from V. 



For any tree t embeddable in G at v, 7r„(t) is a DA on G with source {v}. For any finite DA 
A^'"^ with source {v} in G, we set 

A,(^M):= ^ (9) 

t : 7r„(t)^A{''} 

where the sum is taken on the (necessarily finite) set of trees embeddable in G at v. Our proof of 
Theorem 2.11 consists of proving the two following equalities valid for any finite DA A with source 

D^{A) = A^{A) and A^{A) = Xv{A). (10) 
The key lemma in the proof of the first equality is 

Lemma 2.14 For any finite DA ^{"> with source {v}, we have 

y: (-1)1*1 = (-i)i^^^^i 

t : 7r„(t)=A{''> 

Proof : This is true if A^'"^ = {v}. Assume that this is true for all DA having k cells or less than 
k cells. Take a DA A^""^ with k + 1 cells. Among those cells, there exists w such that w has no 
child in A^''^ Consider the DA A' = A^""^ \ {w} and compare the set of trees ^ whose embedding 
is A^'"^ with 4' the set of trees whose embedding is A'. The trees in ^ are obtained from those of ^' 
by addition of leaves which embedding is w. Hence by the natural projection from ^ into (see 
Fig. 6) 

C = U ^* (11) 

where Kt is the subset of ^ containing the trees obtained from t by the addition of some leaves. 
The Kts form a partition of ^. One has, for any t G 

But given t, there is a maximal tree in Kt (for the inclusion) obtained from t by the addition of 
a set Rt of leaves. The others trees of Kf are formed by t together with a non-empty part of Rt. 
Hence 

^(_1)M= J2 (-1)1*1+1^1 = (-i)l*l((i-i)l«*l-i) = (-1)1*1+1. 

reKt BMcRt,\B\>0 

By (11), X]T-g^(— 1)'^' = — Z^tg^'(— 1)'*' f'^d we conclude by recurrence on the number of cells of A. 

□ 

Now, we conclude the proof of Theorem 2.11. 
Proof of Theorem 2.11. On one hand, by Lemma 2.14, Dv{A) (defined in (5)) satisfies 

D.{A)= E ( E (-1)'*'+M= E (-1)1*1+^ = A.(A). 

B:B^A,S{B)={v} \t : 7r„(t)=S / t : ■Ky{t)-<A 



-t n 




Figure 6: On the first picture, a DA A^"^'^ (whose cells are colored). On the second one, a DA A' 
obtained from A^^^ by the suppression of a cell having no child. The three trees with straight lines 
represent trees embeddable in A' at v. The trees with in addition some of the doted lines are the trees 
embeddable in A^^^ at v. The square nodes are the nodes whose embedding are w 



We now establish that A„ owns the same recursive decomposition (3) as If |^| = 1 then 
Ay (A) = 1 = Xvi^)- Assume that the source v of A has d children vi, . . . ,V(i- For any i £ {1, . . . ,d}, 
we have 

A{A^^):=- (-1)'*'' 

t : 7r„;(t)^A"i 

where A^^ is the maximal DA included in A with source Vi. Decomposing the trees t embeddable 
in A according to their embeddings in the A'"\ we get 

AMI = im>o\ E U^-if 

\Bc{vi,...,va}i(^B 

= vi>o( E (-i)i''in^'''(^''')) = vi>on(i-^'''(^''))' 

\Bc{vi,...,va} i(^B / i=l 

Hence A^, owns the same recursive decomposition (3) as Xv, and finally A^, and Xv coincide on the 
set of finite DA with source v. □ 

2.5 Why is the embedding of trees in DA "natural"? 

According to (1), we have 

Ep{X,)=EplB^{p) J] (l-^c)J (12) 

\ c: children of a; / 

where B^{p) is the Bernoulli (p) random variable llc^=a (recall that all the random variables B^[p) 
are independent). We now expand the right hand side of (12), conditioning by A^^^ : 

Ep{X,)=Ep(ep(b^p) II (1-X,)|aWJJ. 

\ \ c: children of a; / / 



^ A 



Assume now that A^^^ is known (and finite). By construction, B'^ijp) = 1 for any u in A^^^ and 
B^{p) = for any u G 'P(A^^^)- We now replace each of the X(. by a product involving its children 
B^ip) Wf' children of Jy^ ~ ^c')i recursively; we stop this expansion when the children of a node are 
all perimeter sites, since in this case all the X^i equal 0. 

When one expands this formula (there are no more variables X^, since they have been replaced 
by those concerning their children, until the perimeter of A'^^^ has been reached), one builds a 

tree : we have ric'x' children of c(l - ^c) = Ecc{ children of c} ric'GC -^C corresponding in the tree 
like decomposition to the choice of some 1 and some —Xc'. The "1" in the decomposition are the 
leaves of this tree, and when C is chosen, \C\ is the degree of an internal node in this tree. One then 
sees that the final value of Xr^ knowing A^^^, is X]f7r(t)^A{'^}(~-'-)'*'^^ other words A^^A'^^^). 
Indeed, the terms (—1)1*1+^ counts the number of times —X has been chosen in the expansion : 
this number equals the total number of children in the tree, |t| — 1. Hence, for p < Pgrit' 

Ep(X.) = 1:^^(1 -P)l^(^)l (-l)l*l+'=IEp(A.(AW)) 

A t:7r(t)^A 

(x) 

On the other hand for p < Rx < Pcnv 

-G.(-p) = ^p^(l -p)l^(^)l (-1)'^'^' = MDx{A^'^)). 

A A'<A 

Then, -D(^) = A(^) = x(^) is indeed an explication at the level of objects of — Ga:(— p) = Ep(Xa;). 

2.6 Directed animal with compact sources 

The previous subsections deal almost only with DA with sources of cardinality 1. But most of 
the results stated there admit some generalizations to DA with larger sources. The main results 
that we want to state is the following generalization of Theorem 2.6. 

Theorem 2.15 Let G = {V, E) he an agreeable graph, k he a positive integer, S = (si, . . . , s^) he a 
free finite suhset of V and p smaller than the radius of convergence of Gg . Under ¥p, the random 
variables {Xx)xeS o,f^ well defined by (1), and we have 

^p[\{xA =Pp(X, = l,xG5) = (-l)l^lGf(-^)). (13) 

Observe that 

Pp(X, = l,xG5)= ^1-^1(1-^)1^(^)1 (14) 

A,AG.A(5,G),|A|<+oo,Xx(A)=l,a;G5 

for any p G [^■>Pcrit) "-"^ occupied on their source S. 

Proof . The first equality is clear, and to prove the second one a generalization of the proof of 
Theorem 2.11 is needed; we skip the details and give only the main lines of this generalization. 

First, following the lines of Proposition 2.8, we have 

-Gf(-p) = ^,{Ds{A')) (15) 



where for any finite DA A, 

DsiA) = ts^^)=s E (-1)1^'!+^ (16) 

A' : A'~(A,S{A')=S 

Let f = (ti, . . . , tfc) G T*^ be a forest with k trees. We say that f is embeddable in G at if for 
any i e {I, . . . ,k}, tj is embeddable in G at Sj. We set 

^5(f)= U TTsM 

i(^{l,k} 

the union of the embeddings of the tj's on the s'^s. For any forest f embeddable in G at S, TTsi^) 
is a DA on G with source S. For any DA A with source S in G, we set 

f=(ti,...,tfc) : ferfc,7rs(f)-;A 
where the sum is taken on the set of forests embeddable in G, and 

As{A):=ls^^^^S E ^^(^')- (17) 

A' : A'~(A,S{A')=S 

A simple analysis yields Ys{A) := Yl^^i A^. (Aj) where A^ is the maximal sub- DA of A with source 
Sj. For A = S, we have l5(^) = 1. Adapting the proof of Lemma 2.14 to forests, we get 

F5(A) = (-1)I^H5| 

By (16) for any A with source S, 

Ds{A) = (-l)i+l^lAs(^). (18) 
It remains to relate As{A) with HsgS^'^s- definition of the gas model of type 1, we have 

{to, Xs(ijj) = 1, s & S} = {io,CsiLo) = a, s E S, Xy^Lv) = 0,y G Si} 

u,S{A'{co)) = S,l[(l- XyiA^'H^))) = 1 

where Si is those vertices in G, children of elements of S. Write 

{u;,As{A^{u;)) = l} = {a;, C,(a;) = a, s G 5, J] A.,(4^'^(u;)) = 1} 

sies 

Now, consider Si{i) the children of Sj in G. By the proof of Theorem 2.11, 

A,,(Af')(c.)) = l{,,c.,H=a} n a-^yiA'{co))) (19) 

2/eSi(i) 

and then 

n A,,(4^^)(a;))= n (l{a;,a,H=a} n (1 - A,(^^(a;))) 



Now, notice that a node y in Si may be in several S'i(i)'s, but since 1 — /S.y{Ay {u)) G {0, 1} the 
repetition of this factor in the product (19) may be simphfied; we then obtain 

sies \ yeSi J 

Using (10), this is equal to t{u,,Cs(u,)=ayseS} Uy^sS^ - XyiA^'^K^))- Finally we get 

{u;,Xs{u;) = l,seS} = {uj, AgiA'^ (oj)) = 1} 
and this, together with (18) and (15) allow to conclude. □ 



3 Gas model of type 1 and DA on 

3.1 New derivation of the GF of DA counted according to the area 

In this part we work on the square lattice Sq = where as said in the Introduction, Sq is 
viewed as a directed graph where each vertices {x,y) has children {x,y + 1) and {x + l,y + 1). 
We represent this lattice as on Fig. 1. The real number p^^^ • — Pcrit 

does not depend on x by 

symmetry, and for any finite source S, p^^.^ = P^rif According to Proposition 2.2, p^^jj > 1/2 and 
by Proposition 2.4, for p < p^^^f the gas occupation is Pp a.s. defined everywhere. The radius 
can also be shown easily to be positive using a simple surjection from the set of DA with n cells in 
Sq into the set of trees having only internal nodes with total degree 2 and 3 and n nodes. 

The set Lj := {{x,i — x),x G Z} is called the ith horizontal line and we denote by Zj its gas 
occupation : 

Zi{x) = X(^^^i_^), for any a; gZ. 

By construction Zj is a process indexed by Z taking its values in {0, 1}. It is invariant in the strong 
sense : for any A; G Z the process {Zi{x))xez and {Zi{x + k))xez have the same distribution. By 
construction also the processes Zj and Zj+i have the same distribution. Since Zj+i is built from 
Zi via the gas evolution (1), for any i G Z, Z, and Zj+i are related by : 

Zi{k) = Blip) (1 - Zi+i{k)) (1 - Zi+i(fc + 1)) (20) 

where the B^p) are independent and identically distributed (i.i.d.) Bernoulli(p) random variables, 
independent of Zj+i. 

The existence (and the construction) of a solution is guaranteed when p G [0,^^^^^) by Proposi- 
tion 2.4, and a solution, say Z* is characterized by the equation (14) : let S = {{x,—x),x G /} be 
included in Lq, with / finite. We have 

¥{Xs = i,seS) = ¥{z*{l) = i,lei)= Yl pl^l(i-p)l^^^)l. 

A,AeA{S),x4A)=l,seS 

By inclusion-exclusion, this determines the finite dimensional distribution of Z*. 



Proposition 3.1 For all p G [OiP^rit)' ^^^'^s exists a unique law fi of processes taking their values 
in {0, 1} solution of (20) (in other words, there exists a unique law /x such that if the process ^i+i 
is distributed then so do Zi). 

Proof. The uniqueness of the solution wih be proved to be a consequence of Theorem 2.15 : 
assume that Z'q has a certain law /i' which is solution of (20). One then builds an infinite sequence 
(Z-)j<_i such that for any Z[ is obtained from Z[j^^ by 

Z[{k) := Bf{p) (1 - Z[^,{k)) (1 - Z[^,{k + 1)) for i < -1 (21) 

where B^\p),i < —l,k G Z is an array of i.d.d. Bernoulli(p) random variables. 
Consider / a finite subset of Z. We will examine the probability ¥(Z'-{1) = 1, Z S /) and show that 
it converges when i +oo to the distribution of F{Z*{1) = 1,1 £ I). Since P(Z-(/) = 1,1 e I) 
does not depend on i this entails that the finite distribution of Z' and those of Z* are equals. This 
result being valid for any finite set I, this implies that Z' has the same law as Z*. 

Let / be fixed, and consider S{i) = {{x,i — x),x G 1} the subset of the ith line Li, which is 
simply obtained by the translation of S := S{0). Under is a.s. finite; this implies that for 

any e > there exists such that PpdA"^! > n^) < e. 

Consider now S{—n^), the subset of L-n^- We have Pp(| A'^(~"=)l > rig) < e, by invariance by 
translation. When |A'5'(-"-)| < n^, (A^(-"-) U P(A^(-"-'))) n Lq = 0. In this case, the values of 
{Xs)seS{-n^) does not depends on the value of X on Lq (and then neither on its distribution on 
this line), since as said in Section 2.3, the gas occupation on a subset S is a deterministic function 
of A"^ and its perimeter sites (and then does not depend on the other sites). 

The formula (21) defining is the same as that defining the gas model of type 1, and then, 
if A'^*^""''-' satisfies |A'^(~"=)| < n^, then (Z'_^^{l),l £ I) does not depend on Z'q. Therefore, for any 
n > rig, we have 

|Pp(Zi„(/) = 1,1 el)- FpiZ^l) = l,lGl)\<e 
and then Pp(Zi„(/) = l,l£ I)^ ^p{Z*{l) = 1,/ G /). □ 

ZL2 

zU 
zu 

Observed sites does not depend on Zq 

Figure 7: The colored sites are the sites where the Bernoulli(p) coloring = 1. The gas occupation on 
the three first cells on the last row depends only on the DA (with these three cells as over-sources). 
The height of this DA is so small that the gas occupation of these cells does no depend on Z'q. 




Comments 3.2 The sequence (Zi) forms ("a vertical") Markov chain with infinite number of 
states. If one wants to determine the value of Zi{x) using the Bernoulli tossing allowing to re- 
alize this Markov chain, only the DA of calculus A^^''*^^) is needed. This is reminiscent of the 
Propp & Wilson coupling from the past technique to simulate a Markov chain under the stationary 
distribution : only the last transitions of the Markov chain sufficient to determine the current state 
are needed. 

We now determine the distribution ^ of the processes Zi and the density of the gas model of 
type 1. 

Let Z := Zq = (Zo(x))x>ez be the gas process on the hne. We denote by (i?*)j>i (resp. (i?°).j>i) 
the successive sizes of the blocks of consecutive occupied (empty) positions at the right of zero. If 
a block contains some negative positions, only the part at the right of zero is counted as on Fig. 8. 

DO -D9 DO D* DO 

-Dl -Dl -D2 -D2 -D3 

• — o — o — • — • — • — o — • — • — o — o — • 

1 

Figure 8: On this example, the black cells are occupied and corresponds to the places where Z = \, 
and the whites one are empty. 



A random variable G is said to follow the geometric distribution with parameter a, we denote 
G ~ g{a), when P(G = A;) = a(l - af-"" for any k>l. 

Theorem 3.3 Let p G {0,pcrit)- The density of the gas model is given by 

Pp(Z(0) = 1) = ' 
l/ao + 1/a. 

where 



+ Vl + 2p-3p2 -1+P+ ^l + 2p-3p2 

a, = and Oo = , 

2p 1+P+ Vl + %r^3p2 

Under ¥p the random variables B' and B!^ fori > 1 are independent and independent of Z{0), and 
B- r^g{a,) and Bf r^g{ao). 

Hence, for any b°,b' > 1 

P(Zo(0) = X, B° = 6°, 5- = 6-, i G {1, . . . , k}) = n = ^?)^(^' = 

l a, + l/ao 

1=1 

where G' ^ G{a,) and G° ~ G{ao), and a(l) = a, and a(0) = Oq. 

The law fi of Z is characterized by the properties given in this Theorem and the invariance by 
translation, since this characterizes the finite dimensional distributions of /i. We have no combina- 
torial explanation to the properties of the blocks, but these distributions may be guessed from the 
result of Bousquet-Melou cited at the beginning of Section 3.3. 



Comments 3.4 • The proof that we provide below shows that the distribution described in Theorem 
3.3 is a solution of equation (20) not only for p £ {0,Pcrit) but also for any p G (0,1). For 

p € [pcrit, 1) we didn't find an argument to show that there is a unique solution to (20). In the case 
p = 1, the solution is not unique since an alternating sequence of blocks, the occupied block having 
size k and the empty one, k + 1 is conserved by the gas evolution for any k > 1. For p = 0, the 
process Z null a.s. is the unique solution. 

• One may also view Z as built as juxtaposing alternatively occupied and empty blocks with the 
prescribed distributions, starting from — oo. Viewed like this, the density can be computed using a 
simple renewal argument or by the law of large number. 

• Assume that Z\ = {Zi{x))xez and Zq = {ZQ{x))xez have the property announced in the Theorem, 
but that a, and Oo are unknown. To compute a, and Oo we derive and solve the two following 
equations 

(1 — ao)p = 1 — a, and ao = a,(l — ao)p. (22) 

Indeed, a block ofl in Zq is placed "under" a block ofO in Z\. Once a block of 1 is begun in Zq, there 
is an additional 1 at the right of this block with probability (1 — ao)^' (an additional is needed in Z\ 
and the right Bernoulli tossing is needed). This gives the first equation, since in Z\ an additional 
1 occurs with probability 1 — a,. On the other hand, since {x,Zq{x) = 1} = {x,Cx := a, Zi{x) = 
0,Zi{x + l) = O}, the probability ^/^/^'/^^ =¥p{Zo{x) = 1) = p¥{Zi{x) = 0,Zi(x + l) = 0) equals 

to pF{Zi{x + 1) = 0\Zi{x) = 0)F{Zi{x) = 0) = P i/^^^i/a ~ ^o)- ^^^-^ 5^'^^ second equation 
(that can also be obtained by a renewal type argument without using the density). 

3.2 Consequences of Theorem 3.3 

It is worth mentioning that Theorem 3.3 is in fact equivalent to the following representation, 
allowing also to compute the gas density. 

Proposition 3.5 Under Pp, the process {Z{x))xez is a Markov chain with transition matrix 

nZ{l) = l|^(0) = 1) P(Z(1) = 0|Z(0) = 1) \ ^ / 1 - a. a, \ 
1^ P(Z(1) = 1|Z(0) =0) P(Z(1) = 0|Z(0) = 0) J \ ao l-ao) 

under the stationary distribution. 

One may compute from this Proposition and Theorem 2.15 some results about the GF of DA with 
general sources on the square lattice. 

Proposition 3.6 Let S = {si, . . . , s^} where Si = (xj, — a::^) be some points on the principal line 
Lq, such that di := Xj+i — Xi for i £ {2, . . . , fc} are positive integers. The GF of DA on the square 
lattice with source S is given by 



' ' 1=1 



^ ^ a,(l — a, — ao)*^' + ao 
a, + ao 



d2 ds di ds d2 d^ 

Figure 9: On the first picture, the sources considered in Proposition 3.6. The second picture illustrates 
the sources considered in Proposition 3.7. 



Proof. This is a simple consequence of the Markovianity of Z and of Theorem 2.15. Using a 
diagonalization of M (given in (23)), one gets F{Zk = \\Zq = 1) = {M^)^-^ = ■ □ 

For example, if 5„ := {(i, = 1, . . . ,n} then G|^J-p) = ^^-1Z^(1 - a')"-i(-l)-. Then, 
the GF of DA on the square lattice with compact sources satisfies 

y i-v) = , -^^"^ , V(l - a')"-i(-l)" = 

n>l / • I / o J- 

A combinatorial explanation of this formula is given in Gouyou-Beauchamps &: Viennot [11]. One 
may also compute the GF of DA having their sources on different lines. For example : 

Proposition 3.7 Let S = {si, . . . , s^} where Si = (x,, — y?) he some points of Sq such that yi = 
Xi + i (see Figure 9 (ii) ) and di := Xj+i — Xi for i G {2, . . . ,k} are positive integers. The GF of DA 
on the square lattice with source S is given by 

l/a, + tto -^"V a, + Q!o 

1=1 

See Figure 9 for an illustration of the considered sources. 

Proof. We use again the Markovian properties of Z but this time on several lines in the same 
time. We want to compute P(Xs = l,s G S). Notice that the cells of S are not ancestors from 
each others. For any s G S, denote by ci(s) and C2(s) the two children of s in Sq (ci is at the left 
of C2); we have 

P(X, = l,s£S)= F{Xs = 1,X,,(,) = = 0,sGS) 

since if = 1 then a.s. Xci(s) = ^C2(s) — 0- We condition on the gas occupation on si, ci(s2) and 

C2(S2) : 

F{Xs = hsGS) = P(X, = 1,sG5\{si}|X,, =1,X,,(,,)=0,X,,(,,)=0) 

xP(X,, = 1,X,^(,,) = 0,X,,(,,) = 0). 

By Markovianity we have P(X,, = l,X,^^„^) = 0,X,^^,^) = 0) = (Af^0i.2M2,2 and 

P(X, = l,seS\ {s^}\Xs, = 1, = = 0) = P(X, = 1, s G 5 \ = = 0) 

F{Xs = l,s£S\{si}) 



Thus, P(X, = l,seS)= P(X, = l,se S\{si})^(M'^i)i,2. Usingthat (M'=)i,2 = ""^^'Z+Z 
we get the result. □ 

Remark 3.8 Let S = {si, . . . , Sk} be a free subset of Sq. The serie Gs may be computed (theoret- 
ically) using the properties of the gas model, since it suffices to compute P{Xg = 1, s G S). This 
can be done by writing a finite sum involving the occupations of the cells that are above the Sj 's and 
under the line Lj with the largest index such that Lj nS is not empty ( inside the region surrounded 
by the dotted lines in Fig. 10), and using the Markovianity of the gas occupation on Lj. 

r----vOr r;--^-^-r---v 

Figure 10: The points stands for the points of S; the dotted lines surround the (finite) set of cells 
where the summation has to be done. 

The Markovianity of Z implies some (weighted) properties of DA taking into account the occupa- 
tion or not of their sources; here, by weighted we mean that here we deal with probabilities, and then 
each DA A has weight j)l^l(l— p)l^(^)t which depends on its area and perimeter. Let z := (zq, • • • , Zk) 
be fixed in {0, 1}. The probability of the event {Z{0), . . . , Z{k)) = z depends only on the four 
quantities nl'(z) = if{i,i G {0 . . . A; - 1}, (zj, Zj+i) = (1,1)}, nl°{z), n°*{z) and n°^{z), defined 
accordingly. Hence P((^(0), . . . , Z{h)) = z) = P((^(0), . . . , Z(h)) = z') for z = (0, 0, 1, 1, 0, 0) and 
z' = (0, 1, 1, 0, 0, 0). In terms of DA with horizontal sources S = {{i, —i)}i=o...5 = {vi}i=o...5 it leads 
to 

^ pl^l(l-p)l^(^)l= J2 p\^\l-p)\nA)\ 

A,Xv^{A)=z„ie{0,...,5} A,XviiA)=zlie{0,...,5} 

A bijection between DA with the considered occupation of their sources and preserving moreover 
the area and the perimeter would be a direct explanation of this identity, but such a map does not 
exist since the enumeration of these two sets of DA according to the area and the perimeter are 
distinct for DA of area 5 : 

A,\A\<5,Xvi{A)='^i'^^{^'-'^} 

and 

J2 = + (2n« + 10u9)t^ + (1«« + 8uy\ 

^,l^l<5,Xt.i{^)=2-,je{0,...,5} 

3.3 Proof of Theorem 3.3 

One may propose several proofs of Theorem 3.3 : following the proof of Bousquet-Mclou [3], 
one may first work on the cylinder as done in Section 5.2. There the process Z is in some sense 



horizontally Markovian, and conditioned to come back at its starting point. Theorem 3.3 may be 
obtained by a passage to the limit on the length of the cylinder (the limit in distribution here is 
just a convergence of the finite dimensional distribution). 

More directly, another method consists in proving first Proposition 3.5 : three steps are needed. 
At the beginning, suppose that Z\ is a Markov chain with two states and 1. Then using that 
is must be preserved by the gas evolution, there is at most one Markov chain satisfying these 
constraints (easy). Let M be its transition matrix. It remains to show that this Markov chain is 
indeed preserved by the gas evolution. This is done via some linear algebra on M, showing that if 
(Zi(i), . . . , Z\(i + A;)) is a Af-Markov chain, so do [Zq(i), . . . , Zi^(i + A: — 1)) (this is done via some 
quite laborious computation). Proposition 3.1 allows to conclude. 

The proof of Theorem 3.3 that we propose below may be viewed as a translation of what is 
said above. But we think that it is in fact different and is more likely to be generalized to other 
models. Another pleasant reason is the natural appearance of the gas density in the beginning of 
the considerations (see Comment 3.4) : using the other approaches, this is not the case. 

Proof of Theorem 3.3. Let (6°)jG{i,...,fc} ^^^d (&*)jg{i,...,A:} be two fixed vectors in {1, 2, . . . j'^. For 
i in {0, 1} we denote the block sizes of Zj by B° ■ and We assume that .Z^i(O) and the block 
sizes of Z\ have the distribution described in the Theorem, and we prove then that ^o(O) and the 
blocks of Zq have the same distribution. By Proposition 3.1 this is sufficient to prove the Theorem. 

To compute P(Zo(0) = 1, B°q = b°, B'q = b' ,i £ {1, . . . , k}), we will sum on all possible block 
configurations of Zi on the positions {0, . . . , 1 + Yli=o K + We recall that by the gas model 
of type 1, the process Zq is solution of (20) (see Fig. 11). Thus, knowing k consecutive occupied 
cells in Zq allows to guess k Bernoulli random variables and a block oi k + 1 empty positions in Zi 
(that maybe included in a larger empty block). In particular, if Zq{Q) = 1 then ^i(O) = 0. 

The empty blocks of Zq are a little bit more tricky to handle : for each empty block in Zq with 
size 6°, there are h° — 1 cells in Z\ that can be occupied or empty (see Fig. 11). We call these cells 
uncertain cells. 

Notice also that {Zo(0) = 1, i?°o = S'q = h^i G {1, . . . , A;}} = {Zo(0) = 1, i?°o = Bio = 
6*, i G {1, . . . , k}, ^o(X]i=i b* + b°) = 1} since knowing the size of the 2k first blocks and the values 
of Zq(0) implies the occupation of the cell that follows. Hence, we have 

k + 1 




Figure 11: A block of fc occupied positions in Zq determines fc + 1 empty position on Zi (surrounded 
by rectangles). The little squares figure out the uncertain cells, that can not be determined just in view 
of the gas occupation in Zq. 



P(Zo(0) = l,i?°o = ^i,^'o = bl^^ {1, ■■■,k})= , J^"], [{l-a.)pp=^'' (USk) (l-«o)p. 

The first factor i/^^+i/a, right hand side comes from Zi{0) = 0. The other term are 

computed successively, as transition probabihties : 

- the term (1 — ao)p corresponds to the creation of an occupied cell x in Zq knowing that x is 
empty in Zi : 

P(Zo(x) = 1, Zi{x) = 0) P(^o(a;) = 1, ^i(x) = 0, + 1) = 0) 



P(Zo(x) = = 0) 



P(Zi(x) = 0) P(Zi(x) = 0)) 

= P(Zo(x) = l\Zi{x) = 0, Zi{x + 1) = 0)P(Zi(x + 1) = 0\Zi{x + 1) = 0) 

= p{l-ao). 



In other words knowing that x in empty in Zi it is occupied in Zq if there is a transition empty- 
empty in Zi between x and x + 1 and a favorable coloring = a of a; (proba. p). There are 
Ya=i K ^^c'^ terms. 

— The last term (1 — ao)p comes from the creation of an occupied cell after these blocks in Zq. 

- The factor Sij° corresponds to the contributions of the uncertain cells in Zi above the 6? empty 
cells of ith empty considered block of Zq : more precisely 

Si = P(Zo(l) = 0, . . .,Zo{l) = 0,Zi(/ + 1) = 0|Zi(l) = 0). (24) 

Above the empty cells 1, . . . ,1 in Zq there are I — 1 incertains cells in Zi, at position 2,...,l. 

For I = 1 there is no uncertain cell above and we find 5i = (1 — ao)(l — p)- For A; > 2, we 
write Sk = 8^ + S^, where 5* (resp. S^) corresponds to formula (24), where in the RHS is added 
Zi{l) = 1 (resp. Zi{l) = 0), in other words where the occupation of the last uncertain cell is 
specified. For = 1, we take the convention 8° = Si and 5* = 0. For k > 2, there is a simple 
decomposition of S^ and S^ obtained by removing the last cells of this block: 

' S; = S,°.i(l -ao){l-p) + Sr.i(l -£.„)(! -P) 

= ^'-' (l-!:)Tl-p) +^'-'^-''-> 

Using (22), this finally gives 5";° = (1 — p){l — ao)Si-i and S' = (1 — ao)pSi-i, then Si = {1 — 
aoY-^Si = (1 - Qo)'(l -p). Thus 

P(Zo(0) = 1, = bl Bl, = 6-, i G {1, . . . , k}) 

1/ao + l/a, 

using that (1 — ao)(l — q;,)(1 — p) = aoa,, and (1 — ao)p/ao = l/a., we get 

P(Zo(0) = l,5°o = 6°,5-o = 6-,zG{l,...,A;}) = . . , {[FiG* = b'nG^ = b°) 

-L / Cto ~r -L/ (X^ . 

' ' 1=1 

which is the expected result. 

The computation of P(Zo(0) = 0,B°q = b°,B'Q = b',i ^ {1, . . . ,k}) can be performed along the 
same lines, except that here the two cases ^i(O) = 1 and ^i(l) = 1 have to be considered. □ 
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4 Simultaneous construction of DA and gas model of type 2 



We construct in this Section a probability space on G in such a way that the density of the gas 
model equals up to some change of variables the area and perimeter GF for DA. This approach is 
very similar to that of Section 2. Here a probability measure having two parameters is needed. The 
idea is again to transfer a quantity revealing the perimeter size via a DA of calculus. The main 
difference with Section 2 is that here, the revealing quantity is random when it was deterministic in 
Section 2 (it was x(^))- We endow the perimeter sites with i.i.d Bernoulli(p) random variables, and 
we will transfer to the source the minimum of those random variables. Assume that N is random 
and that you know the law of M := minjSi, . . . , -Bat} for . . . , i.i.d. Bernoulli(p) random 
variables, independent of N . It is straightforward that the distribution of M characterizes that of 
A^. This is morally what wc use here. 

Consider G an agreeable graph, and the probability space = {a,b,c]^, endowed with the 
probability measure 

where PaiPbiPc arc three non negative parameters summing to 1. In other words, wc have again 
a random coloring of G, where the nodes colors are i.d.d. and follows the following distribution 
(we write C* the color of x). Under ^pa,pi,pc: fo^ x, is a with probability Pa, C* is b with 
probability pb and C* is c with probability pc- 

We define now Af (w) in the same way as A'^(a;) is defined in the beginning of Section 2.1 : 
Af(a;) has again as set of cells the maximal DA with over-source S whose cells are all a-colored; 
now the perimeter sites of may have the color b or c. In order to take into account these colors, 
we introduce 

Vb{A) = {x,xe V{A), CI = b}, and Vc{A) = {x,x e V{A), C* = c}, 

the perimeter sites of A having color b and c. The set of DA on G with bi-colored perimeter sites 
and source S is denoted by A*{S, G), and wc add in the same way a star to the variables already 
defined in the previous section to make visible that we are working with colored objects. 

The following proposition is the analogous, or more precisely a refinement, of Proposition 2.1 : 

Proposition 4.1 Let G = {V,E) be an agreeable graph. 

(i) Let B G A*{S,G) be a finite DA with source (exactly) the free set S. We have 

'^Pa,Pb,Pc{-^-k —-^J—PaPb Pc 

{ii) Let B G A*{S, G) be a finite DA with over-source the free set S. We have 

'^Pa,Pb,Pc\-^* — ^ I — Pa Pb Pc ) 

where for d G {6, c}, Vs^iB) = VdiA) U{x,xe S, C* = d}. 



Our gas model of type 2 on G is built as follows. For any x eV and u eU, set 



if C^{oj)=b 

1 if C^{u) = c (25) 
^ X^{uj) if C^{uj) = a 

^ c: children oi x 

= lcj=c + lc5=a^ ^.min ^^(a;). (26) 

^ ^ d: children of x 

Once again this recursive definition allows to eventually compute X* for any x, if Pa < Pcru- 
By a simple recursion one easily notices that the values X* all belong to {0, 1}. 

Remark 4.2 Since the values X* all belong to {0, 1}, the min operator coincides with the product, 
and also with the "and" operator, interpreting the gas occupation as Boolean variables. Bousquet- 
Melou [3] when she considers the case {pi,P2,P2,P2) in the square lattice case works under an 
equivalent model even if she uses a different vocabulary (we refer to Section 5.2 for some hints). 

Theorem 4.3 Let G = {V^ E) he an agreeable graph. 
(i) For any x E V and Pa < P^}t we have 

^Pa,V„vSK) =PC + ^Pa,P„pMn{A.*)\ =0)=Pc + G-(p„,Pe), 

where G^{u,v) = \A\>o S(A)=x'^^^^^^'^^'^^^ of the set of DA with source {x} counted 

according to their area and perimeter. 

{a) Let S be a free finite subset of G, and Pa < P'crit- ^Pa,Pb,Pc' we have 

^Pa,P„P. ( n J = ^Pa,P.,PciK = l,XeS)= G'{pa,Pc). (27) 

where G'^{u,v) = XIa <S'(A)c5 of DA with over source S counted 

according to their area and perimeter. 

Proof {i) The first equality follows the fact that X* = 1 if and only if C* = c or A^^-*^ contains x 
and has only c-colored perimeter sites. For the second equality, write 

Pp.P.P.(in(Al^^)| = 0) = Yl pI^^P^^^^^^ = G^iPa^Pc)- 

A,\A\>0,SiA)=x 

The proof of (a) follows the same lines : if all the values X* for x E S equals 1, then the DA ^4"^ 
has all its perimeter sites c-colored, and the points of S* \ S", must be c-colored. □ 

4.1 Comments 

The gas model of type 2 on the square lattice (or on an other lattice), under Fp for p < pcrit is 
well defined, and one may again study the process X^ when x traverses a choosing line. We can 
establish a formula similar to (20) and a Proposition similar to Proposition 3.1 for the gas model 



of type 2, using the same arguments. In order to compute the density of the gas model of type 2 
(this should provide a way to compute Pcr-it which is unknown) a result analogous to Theorem 3.3 
is needed: a description of the gas model on a line, or on "something". But, the gas model of type 
2 even if very similar to that of type 1 has a very different behavior. Except for very special values 
Pa,Pb and pc, it is not Markovian on horizontal line (and no Markovian of order 2, and we think 
no Markovian of order k for any k). 

We want here to point out that some gas models generalizing the gas model of type 2 may have 
the wanted property to be Markovian (or have some suitable structural properties). An idea would 
be to enrich the gas in building a model of gas taking its values in a set larger than {0, 1} (the gas 
model of type 2 will appear as a kind of projection of the generalized model), and again use the 
min operator inside the DA of calculus. The minimum of some random variables with a prescribed 
law f would replace the minimum of Bernoulli random variables. By chance, maybe there exists a 
parameter u (or a family of distribution i/) for which the gas has a simple description. 

5 Appendix 

5.1 Proof of Proposition 2.2 

Consider a random DA A'^. Denote by Lq = S Ci . Then for any i > 1, set 



The random sequence (Lj) gives a decomposition of A"^ into layers : the cells in Li arc the children 
of the nodes in Lj-i that are not children of any cells belonging to Lj, j < i — 2. We want to prove 
that if p < 1/K then Fp a.s., there exists k such that \Lk\ = 0. 

The number Nc of children in ^4 of a given cell c of ^4 has law Binomial(o(c),p), where o(c) is 
the out degree of c in G, that is 



Therefore, given = li for i < k, the random variable has law Binomial(o(C),]9) where 

o(C) is the number of children of the nodes of Ik in G that are not children of any nodes of 
the li,i < k. Hence, given Lj = Zj for i < k, the random variable is smaller for the 

stochastic order than Binomial(iC |/jt|,p) : we say that X is smaller than Y for the stochastic order 
if ¥{X > x) < ¥{Y > x) for any x (we write X <s Y). This is simple property of the binomial 
distribution clear via the representation of binomial random variables as sum of i.i.d. Bernoulli 
random variables. A simple iteration, shows that \Lk\ is smaller for the stochastic order than Zk 
where Zq = \Lq\, and where given Z^^i = z^^i, Z^ is Binomial(i^2;fe_i, distributed. There exists 
a probability space, on which one may construct two sequences (i^) and (Z^) such that 



Li = {v,v E A ,v has a father in Li-i,v ^ Uj<i-iLj}. 




(iLfcl) and (Zi) {Zk) 



and such that a.s. Z^ < Z'j^. On this space, infjfc, = 0} < T := infjfc, = 0}. 



But, Z'^ is simply a Galton- Watson process with offspring distribution Binomial(i^,p) (see 
Athreya & Ney [1]), and then it is classical that such a Galton- Watson process eventually dies out 
(r < -|-cxD a.s.) if the mean of its offspring distribution is smaller than 1 (if \Lq\ is finite, which is 
the case here); here it is the case if p < XjK. □ 

Comments 5.1 If G is a tree in which at depth k, all the nodes have k + 1 children, thenpcrit = 0. 
Indeed, under Pp, the mean number of children of an individual having depth larger than 1/p is 
larger than 1. In the following levels, the number of individuals are larger for the stochastic order 
than a super- critical Galton-Watson process, for which T = +00 with positive probability. 

5.2 Square lattice DA counted according to the area: The historical approach 

The content of this section 5.2 is a reformulation of a part of Bousquet-Melou [3] which is in 
some sense the mathematical justification of the work of Dhar [10] and concerns the enumeration 
of DA on the square lattice using a gas model of type 1 (this terminology is not introduced there) . 

First, instead of working on the square lattice, a cyclic directed square lattice having n cells 
in each row is introduced : the cylinder Cy(n) = {'L/n'L) x Z, for n > 2 which is the lattice Sq 
quotiented by a congruence relation : (a;,y) has children [x mod n,y+l) and (x + l mod n,y-\-\). 

Consider a DA A with source C in Cy(n). Removing the first row of A, the following formula 




Figure 12: A DA on the cylinder. The maximum source has 4 cells, the area is 10. The right most 
and left most vertices are identified. 

is easily checked, 

Gg:(x)=xl^l ^ni^^^ (28) 

DcAf{C) 

where AA(C) = {i,i + 1 mod n,i £ C} is the perimeter sites of C. 

Then, Bousquet-Melou introduced a gas model equivalent to the gas model of type 1. The rows 
of the cylinder are labeled by integers, the i + 1 row being above the ith (see Fig. 12). The gas 
occupation on a row is random vector X = {Xi, . . . , Xn) taking its values in {0, 1}"; its distribution 
is described by Fc = F{Xi = l,i £ C,Xi = 0,i £ "C) or by fc = P(X, = l,i g C) using the 
inclusion-exclusion principle. For any j, the vector X^ giving the distribution of the gas on the jth 
row is given from the X^~^^th by the following stochastic evolution (gas model of type 1): 

Xf = B{{p){l - Xt')il - Xl+\ „)) for any z G Z/nZ,j G Z 

where the (p) are i.i.d. Bernoulli(p) random variables (in other words, a cell i in the jth row 
is occupied with probability p if and only if i and i + 1 are empty in the j + 1th cell, and if a 
Bernoulli (p) tossing is a success, the Bernoulli tossing being independent). 



One gets /^(p) = pl'^IP(X^ = 0,i E A/'(C)) and by the inclusion-exclusion principle, 

/C(P)=P"'' E (29) 



this in terms of F^'^^ reads 



The sequence (X^ ) forms a Markov chain which is aperiodic, irreducible with a finite number of 
states. Hence, there exists a unique invariant distribution F; the corresponding function / satisfies 

/c(p)=pl^l E (30) 

D<zM{C) 

Now, the system (30) is solved explicitly by checking that 

Fd{p) = ^ (-^) (1 - p)l^(^)l where Z„ = Yl ( T^) (^ " ^')'-^^''^' (^1) 

is solution. Notice that for each C, Fq is a rational function of p since it is solution of a linear 
system with polynomial coefficient in p. Observe (28) and (30). If / is a solution of (30) then 

Ggf(-p) = (-l)l^l/c(p) 

is a solution of (28) (in fact there is only one solution up to a multiplicative constant, which is fixed 
taking — J2c ^c^~P^ ~ -'-)• Denote by Gc^ the series of DA with source c, a unique cell. We have 

-G^y(-p)=/c(p)=P(Xi = l), 

the so-called density of the gas model. The series Gc^, we should write G^^"'\ is the series of DA 
in the cylinder. It is clear that its first n coefficients coincide with that of the series of DA on 
the square lattice with source O. Hence in the space of formal series 

Qcy(n) _^ QSq_ 
n— »oo 

The explicit solution of the gas model on Cy(n) is derived as follows : using that \M{D)\ = 
\D\ + \Mr{D)\ where Mr{D) = {i,i eD,i + l^ D}, reads EDc{i,...,n}P'^'(l - P)'^''^'^^' and is 
guessed (or viewed) to have the following form 

n 

Zn= Y JJV'(yi,yi+i) where 

(s/iv,yn)e{o,i}" «=i 

y(0, 0) = V{0, 1) = 1, V{1,0) = p{l — p), V{1, 1) = p and yn+i '■= yi- This corresponds in some 
sense to a "cyclical Markovian model" where a transition ^ 1 or — > is "counted" 1, a transition 
1 ^ is counted p{l — p), a transition 1 — > 1, counted p. In other words, since yn+i = yi, Zn is 



fr(y") where Y 



; this is equal to A" + where Ai > A2 are the two eigenvalues 



1 1 

p{l — p) p 

of y, that are : Ai,2 = ^'^^'^'^^^'^^ . It remains to express the density. In Bousquet-Melou [3] this 
is done by differentiation of (31) (where some mute variables replace {p/{l —p)) and (1 —p))- One 
may also compute this density by computing Wn/Zn where Wn = (y")2,2 which corresponds to the 
contribution of the configurations where the first cell is occupied (in other words, yi = Un+i = !)• 
With very simple linear algebra, we get Wn = (Ai(Ai — 1) + A2(A2 — 1))/(A2 — Ai). Since Ai > A2, 

Wn/Zn ^ (1 - Ai)/(A2 - Ai) (32) 

which is equal to j^/(^ + ^)- The convergence in (32) is considered in the formal sense. 

Comments 5.2 In fact, the solution of (31) is more or less guessed by Dhar [9] who works on the 
whole lattice. He "notices" that the distribution of particles on two consecutive lines coincides with 
the hard particle distribution of activity p /{I— p) on two consecutive lines. On the cylinder, this can 
be interpreted saying that the probability to see j occupied cells on Gi U Gj+i is \{p/{l — p)y where 
X is a constant. Given this, the distribution on a line is then just a marginal of this distribution, 
easy to compute. 

In [3] , Bousquet-Melou generalized the study of gas model allowing some new evolutions between 
lines. With these tools, she is able to give some formal link between density of a gas and GF of DA 
on the square lattice (and also on other lattices) counting DA according to several parameters, the 
area, the perimeter, the right perimeter and the "loops". These links are once again formal. For 
example, she derived the following formula : the area and perimeter GF for one-source DA on the 
cyclic square lattice is 

^ t\^y^(^)\ = l-x- p{pi,P2,P2,P2) (33) 



where pi = 1 — x — t and p2 = 1 — x. The quantity pipi,P2,P3,P4:) is the density of a model of gas 
obtained by the rules evolution given in the following table. 



Occupation of x and x + 1 m Zi 


and 


1 and 


and 1 


1 and 1 


Occupation of x in Zq 


Bernoulli (pi) 


Bernoulli (P2) 


Bernoulli (pa) 


Bernoulli (p4). 
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